Quantum computation of multifractal exponents through the quantum wavelet 

transform 
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We study the use of the quantum wavelet transform to extract efhciently information about the 
multifractal exponents for multifractal quantum states. We show that, combined with quantum 
simulation algorithms, it enables to build quantum algorithms for multifractal exponents with a 
polynomial gain compared to classical simulations. Numerical results indicate that a rough esti- 
mate of fractality could be obtained exponentially fast. Our findings are relevant e.g. for quantum 
simulations of multifractal quantum maps and of the Anderson model at the metal-insulator tran- 
sition. 
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I. INTRODUCTION 

It has been realized over the past twenty years that 
the specific properties of quantum mechanics enable to 
conceive new ways of treating and manipulating infor- 
mation (for a review see e.g. In particular, the 
idea of a quantum computer has been put forth. Such 
a device would perform computation on quantum regis- 
ters usually thought of as made of qubits. It has been 
shown that quantum algorithms can be devised which 
are asymptotically faster than classical algorithms, the 
most famous examples being the Shor algorithm which 
factorizes numbers exponentially faster than any known 
classical algorithm and the Grover algorithm which 
searches a database quadratically faster than any possi- 
ble classical procedure Q ■ 

However, not so many efficient quantum algorithms 
have been found, and it is still unclear which problems 
can be treated faster on a quantum computer and how 
efficiently. One of the first possibilities to be put for- 
ward was the simulation of quantum systems on quan- 
tum computers, which was first envisioned by Feynman 
[1| and then made more precise in several subsequent 
works 0, i, 0, i, i, [13, [IH- Some of these algorithms 
have been experimentally implemented in few-qubit sys- 
tems such as in [l2j . Recent more mathematical works 
have established more rigorously that a quantum com- 
puter can indeed simulate efficiently (i.e. exponentially 
fast) a wide class of quantum systems [H, hj. How- 
ever, these works were mostly concerned with the pos- 
sibility of efficient simulation of a quantum system. To 
be complete and comparable with a classical algorithm, 
a quantum algorithm should not only perform a com- 
putation in an efficient way, but also be able to extract 
information from the result of this computation in an effi- 
cient way. Thus a quantum algorithm should also include 
specification of how information is extracted from the fi- 
nal wavefunction at the end of the simulation. Several 
proposals have been made in order to extract efficiently 
information from a quantum simulation, looking at the 
fidelity [l^, the spectral statistics [3, the localization 



length [TtI. t he Wigner function [T^], or the diffusion 
constants [ll|, [l^ . It has been found that the final gain 
compared to classical simulation depends on the choice of 
the observable and on the measurement procedure, and 
can dramatically change the efficiency of the quantum 
simulation. It is therefore important to explore in more 
detail which quantities can be extracted from a quantum 
simulation and with which efficiency. 

The quantum systems studied up to now in this set- 
ting were in general localized or extended in the com- 
putational basis. However, there exists another class of 
quantum systems intermediate between these two types, 
whose wavefunctions are multifractal. Such properties 
appear in physical systems, for example in wavefunctions 
of electrons in a disordered potential at the Anderson 
transition between metal and insulator [l^, [l^, [2l|, , 
or at the quantum Hall transition [2^. It has been 
shown that simple systems displaying such properties can 
be simulated exponentially fast on a quantum computer 
[13, [13 ■ It is therefore interesting to assess if the mul- 
tifractal properties can be extracted from the resulting 
wavefunction at the end of the quantum simulation, and 
with which efficiency. 

In this paper, we explore different strategies to mea- 
sure the multifractality, and more precisely multifractal 
exponents, from a quantum wavefunction produced by a 
quantum simulation, and assess their efficiency. In par- 
ticular, we assess the interest of the wavelet transform to 
perform such tasks. This transform is a generalization of 
the Fourier transform using basis functions localized in 
both position and momentum instead of the sinusoidal 
waves of the Fourier transform. The wavelet transform 
has been used with great success in data treatment and 
data compression, and has been included in compression 
standards such as MPEG4. It has been shown that a 
quantum wavelet transform (QWT) can be implemented 
efficiently on a quantum computer [1^ [13, [ll], but it 
has been seldom used i n q uantum algorithms, with few 
exceptions 0, [13, [Ml, [U [s^] . As concerns multifrac- 
tal analysis, recent theoretical progress in classical data 
treatment have shown that the wavelet transform can be 
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used as a versatile tool to explore multifractality of var- 
ious phenomena, such as DNA sequences, turbulence or 
cloud structure |34l |. 

Our results show that, although the methods of classi- 
cal multifractal analysis cannot be directly implemented 
efficiently on a quantum computer, suitable modifications 
of the method combined with amplitude amplification en- 
able to get multifractal coefficients from a wavefunction 
with polynomial efficiency. This translates into complete 
algorithms of quantum simulations including measure- 
ments which are polynomially faster than corresponding 
classical algorithms. We also give numerical indications 
that the wavelet transform may enable to assess a rough 
measure of the degree of fractality of a wavefunction with 
exponential efficiency. 



II. MULTIFRACTAL QUANTUM STATES 

Multifractal quantum states are characteristic of cer- 
tain systems intermediate between chaos and integrabil- 
ity. They have been found at the Anderson transition be- 
tween metal and insulator, i.e. between extended states 
and localized states, and also at quantum flail transitions 
[H, [13, im, m, More generally, they correspond to 
a whole class of systems whose spectral statistics, for ex- 
ample, are of the semi-poisson type [s^l intermediate be- 
tween Random Matrix statistics (associated with ergodic 
states) and Poisson statistics (associated with integrable 
or localized states). 

Multifractality properties of wavefunctions are de- 
scribed by a whole set of generalized fractal dimensions 
Dq. For a vector = X]il=i "^iN) in an iV-dimensional 
ffilbert space, the spectrum is defined through the 
scaling of the moments 



N 



(1) 



The multifractal exponents, or generalized multifractal 
dimensions, are related to the spectrum by the relation 



(2) 



In order to investigate properties of multifractal quan- 
tum states, simple models have been devised which ex- 
hibit such properties [2l|, Is^, [sj . A particularly simple 
example is the quantization of a map on the torus which 
has been shown to exhibit a whole range of multifractal 
properties, depending on one parameter. The classical 
map is given by 



p = p + 7 (mod 1) 
q = q + 2p (mod f ) 



(3) 



The quantization of this map yields a unitary evolution 
operator acting on a Hilbert space of dimension N — 



f/(27rfi) which can be expressed in momentum space by 
the N X N matrix (sl, [331 
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with 



-2'Kp'^ /N . From this quantized map one can 



construct an ensemble of random matrices, taking (j)p as 
independent random variables uniformly distributed in 
[0,2^[[33. 

This map has different properties depending on the 
parameter 7. fndeed, for irrational 7, the map possesses 
the characteristics of systems displaying quantum chaos, 
with eigenstates ergodic over the phase space and spec- 
tral statistics of the eigenphases of the evolution operator 
following random matrix predictions. In contrast, for ra- 
tional 7 = ni/n2, the spectral statistics are of the semi- 
Poisson type intermediate between those of integrable 
and chaotic systems (37| . The eigenstates in momen- 
tum representation display multifractal properties stud- 
ied in [38[ . The fractality depends on rt2 and is stronger 
for smaller 712. Thus the set of quantum maps U with 
rational 7 gives a random matrix ensemble with inter- 
mediate statistics (ISRM) whose multifractal properties 
are controlled by the parameter 712. Such a map can be 
implemented efficiently on a quantum computer |25j . as 
an exponentially large vector can be evolved through (j4]) 
with a polynomial number of gates. This map will be 
used as a benchmark for numerical simulations in this 
paper, since it represents a simple but non-trivial exam- 
ple of system with multifractal properties depending on 
one parameter and which can be simulated efficiently on 
a quantum computer. 

On the other hand, this system is complicated enough 
to be hard to tackle analytically. In order to test the 
accuracy of our approach for the estimation of multifrac- 
tal dimensions on a simpler system where all exponents 
can be analytically computed, we consider the situation 
where the coefficients of the quantum state are given by 
a multifractal cascade. Multifractal cascades are exam- 
ples of multifractal measures for which generalized di- 
mensions can be calculated analytically. The simplest 
example was proposed in [s^; it is a special case of a 
Bernoulli measure on a two-scale Cantor set [l^]. It can 
be constructed by the following process: one breaks an 
initial interval into two equal parts, attributes a weight 
Pi to one half and p2 = 1 — pi to the other, and repeats 
the process (keeping the same weights {pi,p2) constant) 
on the newly constructed intervals. After k steps, there 
are 2^ intervals and the weight of each interval is of the 
form p\p\~^ for some r, < r < fc. The generalized 
dimensions can be shown to be given by (40j 
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l0g2(pl+Pl)- 



(5) 



In this way one can construct quantum "cascade" states 
of size 2^ whose amplitudes squared are given by the 
weights of the cascade. 
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III. WAVELET TRANSFORMS AND 
MULTIFRACTAL PROPERTIES 

An important procedure has been developed in recent 
years to extract multifractal exponents from a distribu- 
tion. It uses the wavelet transform [4]| . a generalization 
of the Fourier transform which expands a function on the 
wavelet basis instead of the Fourier basis. Contrary to the 
sinusoidal waves which compose the Fourier basis, which 
have specified frequencies but are extended in position 
space, the wavelets are localized both in momentum and 
position. They can thus probe many properties which 
are difficult to reach with the standard Fourier trans- 
form, such as singularities of the distribution. This has 
made wavelet transforms a popular tool in recent devel- 
opments of e.g. image or sound treatment and compres- 
sion in classical information, such as the formats JPEG 
and MPEG. 

Wavelet transforms are based on a single function g, 
called the analyzing wavelet or mother wavelet. The 
wavelet transform of a function / is a function Tf of two 
variables defined as 



Tf{a,b) = ^ / dxf{x)g ' ^ 



(6) 



Variable a corresponds to the scale at which the function 
/ is analyzed, while & is a space variable. Thus Tf(a,b) 
is a measure of how close the function / is to the mother 
wavelet at point b and at scale a. 

If the function / is sampled as a A'^-dimensional vector 
where N ~ 2", the wavelet transform can be discretized 
and implemented as a unitary transformation, resulting 
in a Fast Wavelet Transform (FWT) analogous to the 
Fast Fourier Transform. The scale parameter a takes 
values 1, 1/2, 1/4, ...1/2"^^, while the space parameter b 
varies over £(a) = {1, 2, . . . , 1/a} at scale a. A discrete 
version of the mother wavelet is constructed recursively 
at each scale. Commonly used mother wavelets for the 
FWT include the Haar wavelet [i^] and the Daubechies 
wavelet [1^. Throughout the paper, we will use the 
Daubechies 4 discrete wavelet transform in the numerical 
simulations. It has been shown that a quantum wavelet 
transform (QWT) implementing such discrete wavelet 
transforms can be performed efficiently on a quantum 
computer [2^, [2^ , namely an exponentially large vec- 
tor X]i=o^V'iN) of size TV = 2" can be transformed in 
a number of operations polynomial in n into a vector 
X^a 6 ^'/'('^j where the wavelet transform at scale 

a is stored on the computational basis vectors |a, .). 

The wavelet transform has been put forward recently 
as a tool for extracting the value of the exponents from 
a multifractal distribution. These exponents are usually 
quite hard to extract numerically, and are very unstable, 
since the values obtained depend on the chosen numeri- 
cal method up to very large system sizes. An example of 
a multifractal distribution (for an eigenvector of ^) is 
shown on FiglU together with its wavelet transform. One 
sees that the wavelet transform does not look especially 
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FIG. 1: (Color online) Example of one eigenfunction of Q 
(top) and its wavelet transform (bottom). The red/grey 
dashed vertical lines on the bottom panel separate the dif- 
ferent scales in the wavelet transform. 



simpler than the original distribution. However, recent 
works have shown [34 . H . |45| that the wavelet transform 
allows to extract the exponents of the distribution /, us- 
ing the maxima of the wavelet transform at each scale. 
Such methods based on the wavelet transform have en- 
abled to extract multifractal exponents in complicated 
physical systems of great fundamental and technological 
importance such as e.g. DNA sequences [i^, fully devel- 
oped turbulence [47| or high-resolution satellite images 
of cloud structure [48j . 

However, these methods use the continuous wavelet 
transform, which is delicate to implement for many sys- 
tems of interest. A variation of the wavelet method de- 
veloped in [i^, IH^I uses instead of the maxima of the 
continuous wavelet transform the sum of the values of 
the discrete wavelet transform, properly normalized at 
each scale. One defines the partition function 



Z{a,q) 



E 

beC{a) 



\Tfia,b)\ 



E 



\Tf{a,b)\ 



(7) 



where a is the scale and C{a) is the interval corresponding 
to each scale a. The asymptotic behavior of the partition 
function at small scales is governed by the generalized 
dimensions as 



(8) 



with Tq — Dq{q — 1). 



IV. QUANTUM ALGORITHMS FOR 
MULTIFRACTAL EXPONENTS 

There does not seem to be a simple way to implement 
efficiently the maxima method on a quantum computer. 
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However, as we show below, the scahng ([8]) of the par- 
tition function ([7]) can be evaluated on a quantum com- 
puter, and gives rise to quantum algorithms for estimat- 
ing multifractal exponents Tq . Our goal in this section is 
to construct quantum algorithms which use Eqs. ©-([l]) 
in order to efficiently extract information about multi- 
fractal properties, in the situation where quantum states 
are obtained by quantum simulation. 



A. Wavelet transform of 

For a given state ip, the weights define 
a normalized measure on the unit interval. Following 
Eq. ([7]), the partition function that describes the multi- 
fractal properties of "0 reads 



Z|^|2(a, g) 



E 

be£(a) 



YjbeC(a) 1^1-0 



^{a,b)\ 



(9) 



where as before a is the scale and C{a) is the interval cor- 
responding to each scale a. The denominator in Eq. ^ 
is a normalization that ensures that at each scale the sum 
yields a proper probability measure. The multifractal ex- 
ponents Tq can be extracted in the limit of small scales 
through the scaling ([5|). As is shown in Fig[2]for cigen- 
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FIG. 2: (Color online) Average of logj ^|^|2 {a, q) as a function 
of the scale a for q = 2 (□), q — 3 (o), q — 4 (A). The slope 
of the fitting straight lines (dashed) gives —Tq. Averaging is 
done over ~ 25000 eigenvectors of Q of size N — 2^* for 
na = 3 (7 = 1/3). 



vectors of (g)), the partition function ^ indeed scales as 
predicted for multifractal wavefunctions, and gives the 
multifractal exponent Tq, whence the fractal dimensions. 
For instance the linear fit for q = 2 has a slope giving 
T2 = D2 = 0.52, which is consistent with the fractal di- 
mension obtained in (38| . 

We now show that it is possible to implement such 
sums on a quantum computer in the case 9 = 2, allowing 
to obtain the multifractal exponent T2 for wavefunctions 
produced by a quantum simulation algorithm. The first 



step is to adapt such a quantum simulation algorithm to 
our present purpose. Usually the vector |-0) we are in- 
terested in comes out of the quantum simulation with its 
components "0^ as amplitudes of the basis states. How- 
ever in ([9]) the wavelet transform is applied to a vector 
whose components are Thus one first has to modify 

the quantum simulation in order to produce j-^ipji) 
rather than X^iV'iN)- To do so, one should build the 
product of two iterates on separate registers, in order 

to obtain a state X^il^o^ S^o^ ^^'^ then use 
amplitude amplification to keep only the diagonal terms 
i = j- Let us detail the steps. We start with two copies of 
an initial state |V'o) = X^iLo^ i^i momentum rep- 

resentation: iV'o) ® IV'o) E«=o^E^"^)^ 

Here ^pQ is the complex conjugate of ipo ^-nd \pi) are the 
basis vectors in the momentum representation. This step 
requires 2?! qubits to hold the values of the wavefunc- 
tion on a A^-dimensional Hilbert space, where TV = 2". 
Let us consider the iterations of such a vector through a 
quantum map. We can apply the algorithm implement- 
ing the evolution operator of the intermediate map U as 
explained in [2^ to each subsystem independently. The 
operator on one subsystem can be described as the prod- 
uct of diagonal operators followed by Fourier transforms. 
On the quantum registers this corresponds to multipli- 
cation by phases followed by quantum Fourier transform 
(QFT). The multiplication by phases of each coefficient 
keeps the separability of the state into its two subsys- 
tems. The QFT mixes only states with the same value of 
the other register attached, and therefore also keeps the 
factorizcd form. Let us see how this works for one iter- 
ation: multiplication of the first register by the diagonal 
operator e"*'^''^ Sij performs the transformation 

i=0 j=0 i=0 j=0 

(10) 

After QFT with respect to pi followed by multiplication 
by e^"'''' the state can be put under the form 

N-l N-l 

E(E^'i*))^f (11) 

Under QFT with respect to q.^ we get 

N-l N-l N-l N-l 

J2iT.^''\p^))^T\pj) - (Ev'.HK»«(Ev^rb^)) 

j=0 i=0 i=0 j=0 

= UH'o)®\ro)- (12) 

Thus iterations can be carried independently on each reg- 
ister. By applying the same steps 2t times, we obtain the 
state ?7*|Vo) ® U*'\ro) = Eto' Ef="o' ® V-f b,)- 

This can be done in 0{tv?) gates if we use the algo- 
rithm of [2^ to implement U . Once this is done, the 
state after t iterations is a product of two copies of 
the state whose fractal dimensions we are looking 
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for. In the momentum representation it takes the form 
j ''Pi'^jlPi) \Pj)- We wish to select in this double 
sum the terms with i = j. This can be done through 
amplitude amplification [5l| , which is a generalization of 
Grover's algorithm The latter starts from an equal 
superposition of N states, and in ^/N operations brings 
the amplitude of a specific state close to one. Ampli- 
tude amplification increases the amplitude of a whole 
subspace. If P is a projector onto this subspace, and 
V is the operator taking |0) to a state having a nonzero 
projection on the desired subspace, repeated iterations of 
V{I-2\0){0\)V-^{I-2P) onV\0) will increase the pro- 
jection. Indeed, if one writes y|0) = PV\0) + {I-P)V\0), 
the result of one iteration is to rotate the state toward 
PV\Q) staying in the subspace spanned by PV\0) and 
{I-P)V\0). Iix= |PT>|0)p, one can check that after one 
iteration the state is {4:X^-3)PV\0) + {4:X^-1){I-P)V\0) , 
with a component along (/ — P)V\0) decreased by Ax^. 

If V is chosen to be [/* (8) [/** and P to be a projec- 
tor on the space corresponding to i = j, the process of 
amplitude amplification from the initial state {ij^o) ® |'0o) 
will bring the probability oi i = j to one. The num- 
ber of iterations depends on the probability inside this 
space compared to total probability. As we want to se- 
lect terms with pi = pj the subspace we are searching 

for has a relative weight X^iLo^ IV'jI^- For a fractal state 
this quantity scales as N^", with a < 1. Thus the total 
number of Grover iterations will be of the order N°'/^. 
So in the worst case of a w 1, the total cost of building 

the state X]fc=o^ IV'iPbj)bj) is ~ t\/N. The more fractal 
the system is, the smaller a is, and the more efficient the 
algorithm is. 

Once the state J2iLo^ \4'i\'^\Pi)\Pi) is built, the second 
register is set back to |0) and the QWT can be applied at 
a logarithmic cost on the first register. This yields a state 
12a b '^W^ ("^j Then measurement of a will give an 

histogram of J^t l^^l^p i'^^ fo^' the different values of a. 
The denominator of ^ can be obtained through classical 
approximations with lower resolution (e.g. half the num- 
ber of qubits), and the exponent can be obtained from 
the slope of logj Z\^\2{a, q). The total cost for applying 
t iterations of the map and extracting the fractal dimen- 
sion is thus 0{tN°'^^) operations, as opposed to tN for 
the classical algorithm. 

Here we considered vectors obtained by iterations of a 
map. It is however important to study also eigenvectors. 
The procedure above cannot be generalized directly, since 
eigenvectors are to be selected by phase estimation first 
before the partition function can be calculated. 

In order to obtain eigenvectors of a map U, one 
first builds the state 2"" X)* l^>|C^*V'o) Et K'>|C^**V5) 
with < i < — 1. The initial state |'0o) should 
be simple enough to be built efficiently (for exam- 
ple |00...0)). To obtain the double sum one starts 
from 2-"J2t\*)\'^o)J2'tW)\'^o)^ easy to obtain from 
10)1-00) |0) l^o) by Hadamard gates. Then one applies U or 
U* in the same way as above on each register t or t' times. 



by using conditional gates, in a manner similar to the one 
explained in [s^l • We then do the QFT on both registers 

\t) and 10, giving 2"" ael^)!^^) a^'l - ^OK'e*'), 
where 9 are the eigenphases of the unitary quantum map. 
This corresponds to performing two phase estimations in 
parallel on each copy of The functions ag are peaked 
around the eigenvalues of U, with \ipg) the corresponding 
eigenvector. Then we use amplitude amplification to se- 
lect the same eigenvalue 9 = —9'. This costs in general at 
most N^fN operations in total. It is possible to improve 
this bound by taking advantage of the fact that the wave- 
function is multifractal, as was the case for iterations. 
For example in the simulation of (|4]) , the wavefunction is 
multifractal in the p representation. So if one starts with 
a basis vector in p representation, the eigenvectors will 
have a multifractal distribution on this particular basis 
vector. Thus selecting 9 ~ ~9' costs a factor A^"/^. The 
selection of the diagonal part of the product of the two 
same eigenvectors costs another N"/^. This implies that 
the total cost to get the same eigenvalue will be A^+", 
as opposed to A^^ for the classical algorithm. 

Then the wavelet transform can be applied at a loga- 
rithmic cost. This yields the numerator of the formula Q 
in A^i+" operations. The denominator can be obtained 
through classical approximations with lower resolution, 
giving a total cost of N^^" as opposed to A^^ for the clas- 
sical algorithm, but with the approximation that the de- 
nominator corresponds to a lower resolution. Then mea- 
surement of a will give an histogram of Eb |2]*|2(a, 6)p 
for the different values of a, from which the exponent can 
be extracted. 

The method described here enables to obtain the nu- 
merator in Q more efficiently on a quantum computer 
than on a classical device. Nevertheless, this procedure 
needs to use a low-resolution approximation of the de- 
nominator in ([9]) which will be obtained on a classical de- 
vice. We have made numerical experiments to see if com- 
puting classically the denominator with twice less qubits 
(thus keeping the polynomial gain for the whole proce- 
dure) leads to a good approximation of the Tg. However, 
we were not able to reach a large enough number of qubits 
to have sensible data with half the qubits, and thus the 
evidence was inconclusive (data not shown), although we 
believe the procedure should give a better approxima- 
tion to Tq than with a computation of both numerator 
and denominators at low resolution. To circumvent this 
problem, it is worthwile to investigate if a modification 
of the partition function ([9|) would be more adapted to 
quantum computation. 



B. Wavelet transform of 

In the preceding algorithm, one had to approximate 
the denominator in ([9]) through low resolution classical 
computation. This entails an approximation of the mul- 
tifractal exponent which can be difficult to control. In 
order to circumvent this difficulty, we studied the ac- 
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curacy of a different formula for the partition function 
which is more efficient to implement on a quantum com- 
puter. This is a generalization of the partition function 
© which involves the wavelet transform of ip instead of 
IV'P, namely 



fc6£(a) 



|T^(a,6)p 



E 



(13) 



The advantage of such a choice for the partition function 
is that the QWT can now be performed directly on a 
state Ei^o^ "01 1*) J without using amplitude amplification 
to build the IV'iP on the registers. This new partition 
function actually also has a scaling a'^« for a — > as 
in as illustrated by FiglJl 
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FIG. 3: (Color online) Average of logj Z^{a, q) as a function 
of the scale a for q = 2 (□), q — 3 (o), q — 4 (A). The slope 
of the fitting straight lines (dashed) gives —Tq. Averaging is 
done over ~ 25000 eigenvectors of Q of size A*' = 2^* for 
na = 3 (7 = 1/3). 



In fact, the quantities J2bec(a) appearing 
in the numerator of Eq. ([TS]) should themselves contain 
some information about fractality of the state ip. Let us 
denote the corresponding multifractal exponents, that 
is the exponents extracted from the partition function 
()13|1 without the normalization, that is from the scaling 



2q 



(14) 



beC{a) 



It might seem at first sight that the computation of t[ 
through QWT could yield an exponential gain over clas- 
sical computation, provided exponentially fast quantum 
simulation of the state is possible. Indeed, for iterations 
of simple vectors (such as basis vectors) by an efficiently 
simulable quantum map such as ([4]), the whole pro- 
cess of iterations of the map and QWT is exponentially 
fast, and the asymptotic behaviour of J2b£C{a) 
can be estimated also exponentially fast by measuring 
a and making an histogram of the probabilities to ex- 
tract the exponent t{. For the multifractal cascade, the 



slopes seem to converge to a precise value which depends 
on the fractality (see Figll]). However, our numerical 
data show that for the quantum map ^ the asymp- 
totic behaviour of Ebe£(a) is independent of 
the system and seems universal. This can be seen in 
FigO where the multifractal exponents extracted from 
J2bec(a) I^i/jI'^j shown. It therefore seems that 
J2bec(a) which can be obtained exponentially 
fast, does not yield useful information in the case corre- 
sponding to actual quantum simulation. We shall return 
to this quantity in Subsection C. 
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FIG. 4: (Color online) t[ as a function of p for the multifractal 
cascade; vectors have size N — 2" with n varying from top to 
bottom from n = 6 to n = 20. 



-0.9 



-1 



-1.1 



-1.2 




FIG. 5: Tq for eigenvectors of ([4]) as a function of q for n2 — 3 
(solid), 722 = 5 (dashed), n2 — 7 (dotted), n2 — 11 (small- 
dot), and 712 = 13 (dash-dot) (ni = 1). Average is done over 
~ 25000 eigenvectors of size 2^^. 



In the remaining part of this subsection, we study the 
quantum simulation of the quantity Z^(a,q), again for 
q = 2. We show that it can be obtained with poly- 
nomial gain over classical computation. Since the de- 
nominator of p3p is obtained exponentially fast, being 
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just J2bGC{a) exponent T2 is now obtained 

without resorting to approximations. 
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FIG. 8; T2 for iterates of column vectors corresponding to 
time t = 1000 of the intermediate map Q computed with 



FIG. 6: (Color online) ra for the multifractal cascade (a, g) ® (hUed symbols) and (a, g) ([HI (empty sym 



computed with Z^^^2{a,q) ((9} (filled circles) and Z^{a,q) 
(|13|) (emptv circles). The values of p are from left to right and 
top to bottom, p — 0.1,0.2,0.3,0.4. The straight horizontal 
lines are the theoretical values T2 = — log2(p^ + (1 — p)^). 



bols) for n2 — 3 (squares) and n2 — 5 (circles) (ni = 1). 
number of vectors averaged for each case is ~ 25 000. 
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FIG. 7; (Color online) T2 for eigenvectors of the intermedi- 
ate map (|4|) computed with Z^^^2{a,q) ([9]) (filled symbols) 
and Z^{a,q) p3|) (empty symbols) for 712 = 3 (squares) and 
n2 — 5 (circles) (ni = 1). Average is done over ~ 25 000 
eigenvectors. The green/grey area shows the corresponding 
range of values obtained in |38i] (top n2 = 5, down n2 =3). 



Figs l6l7l51 show the behaviour of the numerical T2 ob- 
tained from (jl3p for increasing system sizes, plotted to- 
gether with the one obtained from One sees that 
both converge to the same value for large system size. 
In the case of the multifractal cascade (FiglH]), they both 
converge to the analytical value given by Eq. while 
for eigenvectors of the intermediate quantum map (jH) 
(Figl?]) the asymptotic value is within the uncertainty of 
the different numerical methods used in [ssj . For iterates 



of column vectors of ^ (FigH]), one also obtains simi- 
lar results with both methods. This seems to indicate 
that this new partition function gives asymptotically the 
same result as ([9]), and can be used reliably to obtain this 
multifractal exponent. 

To get the quantities X]bG£(o) numer- 
ator of (fTS]) for Z^(a,2), one needs to build the state 
J |T'^,(a, 5)p|a, 5). This can be done through a pro- 
cedure similar to the one exposed in Subsection A. One 
starts from two iterations of the map, followed by two 
QWT. This yields the state 

Y,T^{a,b)\a,b)Y,Ti.{a\b')\a',b'). 

a^b a' .h' 

Then amplitude amplification can be used to select the 
part a = a' and b = b' in this expression, leading in at 
most ViV operations to the state h I^V'C'*' 
From this state, measurement of register \a) yields a value 
of a with probability X]tig£(a) I'^i'i'^^ histogram of 

these probabilities thus yields the exponent. In the case 
of the iterations of vectors through a quantum map such 
as the total cost of the method is at most of order 
t\/N operations, to be compared to tN for iterating the 
vector classically up to time t. Actually the scaling of 
the second moment of the wavelet transform vector 
has an exponent (3 smaller than 1 which means that Ty, is 
also fractal. This gives an actual quantum cost of tN^^"^, 
where (3 < 1 for fractal wavelet transform. 

For eigenvectors of a quantum map, the same reasoning 
as above leads to a quantum algorithm costing of the 
order N'^+°'/'^+l^/'^ operations, as opposed to A^^ for the 
classical algorithm, where a < 1 for fractal wavefunctions 
and (3 < 1 for fractal wavelet transform. 
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C. Possibility of exponential gain 

We have seen at the beginning of Subsection B that 
J2bec{a) can be estimated exponentially fast 

on a quantum computer. The corresponding exponent 
gave useful information for the cascade (see Fig|4]) , but in 
the case of the quantum map ^ , our simulations showed 
that the numerically extracted exponent converges for 
large system size to the same value, whatever the frac- 
tality of the system (see Figl5|). Nevertheless, we present 
numerical data on Figl9] showing that for the multifrac- 
tal quantum map although all exponents converge 
to the same value, the rate of convergence seems to de- 
pend on the fractality of the system. If these numerical 
indications correspond to a generic phenomenon for mul- 
tifractal quantum systems, then a rough estimate of the 
degree of fractality of the system could be obtained ex- 
ponentially fast compared to classical algorithms. 




4 6 8 10 12 14 16 

n2 

FIG. 9: t[ = X^i, l-^V-P foi' eigenvectors of the intermediate 
map (|4]) as a function of n2 (with ni = 1) for different values 
of n (number of qubits) bottom to top n = 6, 7, 8, 9, 10, 11, 12. 
The number of vectors averaged for each case is 25 000. 



V. CONCLUSION 



In this paper we have investigated the use of the quan- 
tum wavelet transform to extract information (namely 
the multifractal exponents) from a quantum wavefunc- 
tion produced by an efhcient quantum simulation on a 
quantum computer. We have shown that various parti- 
tion functions which can be extracted from a wavcfunc- 
tion indeed enable to obtain the multifractal exponents. 
Unfortunately, the partition function which can be ob- 
tained exponentially fast does not yield useful results, 
although our numerical results indicate that it can give 
exponentially fast a rough idea of the degree of fractality 
in a system. We have shown that other partition func- 
tions can be extracted, with a polynomial gain compared 
to classical computation. 

Our results indicate that the quantum wavelet trans- 
form can be applied efficiently to extract useful informa- 
tion for certain types of quantum simulation. They show 
that once the full process of simulation and measurement 
is taken into account, quantum simulations of multifrac- 
tal quantum systems yield polynomially efficient quan- 
tum algorithms, which may be exponential in certain 
cases. Our results also show that the quantum wavelet 
transform is an effective tool which acts in a complemen- 
tary way from the more usual quantum Fourier trans- 
form, enabling to extract certain types of information 
efficiently. 

We thank Marcello Terraneo and John Martin for help- 
ful discussions, CalMiP for access to their supercomput- 
ers, and the French ANR (project INFOSYSQQ) and 
the IST-FET program of the EC (project EUROSQIP) 
for funding. 
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